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FOREWORD 

This report summarizes all work conducted by McDonnell Douglas Astronautics 
Company-East (MDAC-E) in St. Louis, Missouri for the Aero- Astronautics Laboratory 
of the Marshall Space Flight Center (MSFC) under contract NAS 8-26699, "Study of 
Uncertainties of Predicting Space Shuttle Thermal Environment". The Contracting 
Officer Representative for this study was Mr. Homer Wilson and the Study Manager 
was Mr. Robert Masek. Significant contributions to this study were provided by 
Mr. Alan Forney of the NASA MSFC, Dr. Thomas Kane and Mr. Lambert Ebbesmeyer of 
MDAC-E. The transition correlations presented in Figures 48 through 50 were 
generated as a part of independent research and development (IRAD) at MDAC-E 
and are reported herein for information only. 

The International System of Units is used as the primary system for all 
results reported herein. The results have also been reported in the British 
Engineering System of Units which was used as the primary system for all 
calculations made during the course of this study. 
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I. INTRODUCTION AND SUMMARY 

This study was conducted to develop quantitative estimates of the uncertainty 
in predicting aerodynamic heating rates for a fully reusable Space Shuttle system, 
and the impact of these uncertainties on Thermal Protection System (TPS) weight. 
Widely differing temperature predictions for the Shuttle on the part of the several 
Phase A and Phase B contractors showed that a study of this type was required, es- 
pecially when the effect of vehicle weight increases on totally reuseable system 
cost (e.g., $28,000 per pound for the orbiter based on MDC Phase B cost studies) 
is considered. 

The study approach consisted of statistical evaluations of the scatter of 
heating data on Shuttle configurations about state-of-the-art (e.g., Phase B) 
heating prediction methods to define the uncertainty in these heating predictions. 
The uncertainties were then applied as heating rate increments to the nominal 
predicted heating rate to define the uncertainty in TPS weight. Separate evalua- 
tions were made for the booster and orbiter, for trajectories which included boost 
through reentry and touchdown. 

For purposes of analysis, the vehicle configuration is divided into areas in 
which a given prediction method is expected to apply (e.g., lower surface center- 
line, shielded regions, interference regions), and separate uncertainty factors 
and corresponding uncertainty in TPS weight derived for each area. 

Various prediction methods were investigated for application to local areas 
on the Shuttle configuration. These included the prediction methods recommended 
by the Aerodynamic Heating Panel, and alternate state-of-the-art methods. These 
included: (1) direct correlation of the wind tunnel data for a specific configura- 
tion in terms of h/h i.e., normalization of heating data to Fay-Riddell 

ret 

theory, (2) experimental determination of three-dimensional cross flow correction 
factors for candidate theories using data on a specified configuration. 

"Best fits” were made to available wind tunnel data to define the nominal 
heating. The data are assumed to be normally distributed about the "best fit” 
method and multiplication factors corresponding to 1, 2, and 3 standard deviations 
of the data about the fit determined. 

Heat transfer analyses were subsequently made for a completely reusable 
Shuttle system consisting of an aluminum heat sink booster and orbiter with 
reusable surface insulation (RSI) for thermal protection. These analyses showed 
that the largest contribution to the uncertainty in TPS weight for the orbiter 
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occurred on lower surface areas due to the combined effect of heating and boundary 
layer transition uncertainties. The critical regions on the booster occurred on 
lower surface regions and upper nose region, but the major region which consisted 
of LOX and Ll^ tanks were sized by structural requirements. Also, it was found 
that the major contribution to weight uncertainty were uncertainties in wing 
heating because of the large area involved. The booster, on the other hand, was 
much less sensitive to uncertainties in boundary layer transition. 

Although the heat sink booster skin thickness was found to be more sensitive 
to heating uncertainties than the RSI concept, a large portion of the booster is 
cryogenic tankage and these sections of the vehicle are sized by structural rather 
than thermodynamic requirements. Therefore, the weight fraction affected by the 
heating uncertainties is smaller than estimates based on heat sink requirements 
alone . 
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II THE AEROTHERMQDYNAMIC PHENOMENA 

For purposes of this study, the Shuttle mission and configuration concepts 
were assumed to be the fully reusable system developed in Phase B* Figure 1 shows 
in pictorial form the important mission phases. The corresponding trajectory 
variables are contained in Figures 2 and 3. During liftoff and boost the orbiter 
is mated to the top of the booster. While in this configuration, prior to staging 
at approximately 2.135 km/sec (7000 ft/sec), critical booster heating is 
encountered on the forward portion of the top of the booster fuselage due to bow 
shock intersections and channeled flow between the two vehicles. The flow is 
turbulent and (since radiation cooling of the surfaces is difficult due to the 
shading by the other vehicle) temperatures may approach the total temperature of air. 

During the separation maneuver » the aft upper fuselage encounters plume 
impingement heating from the orbiter engines. The booster is powered down and the 
attitude pitched to tt/3 radians (60 degrees) angle of attack and booster reentry 
begun. Booster reentry is accomplished at this fixed angle of attack throughout 
the hypersonic flight regime during which lower surface heating is critical. Due 
to the presence of wings and canards » heating because of shock impingement produces 
local hot spots on the side of the fuselage. Local heating increase may also 
occur on the wings due to nose/wing shock interactions. 

After vehicle separation, the orbiter engines deliver the vehicle to the 
planned orbit, the experiments or other orbital mission is accomplished , and the 
vehicle prepared for reentry. Except for abort or once around missions, the 
boost phase heating to the orbiter is dissipated and the initial TPS temperatures 
determined by the orbit, vehicle inclination , and thermal control provisions. Thus, 
the orbiter TPS is sized primarily by reentry with the possible exception of 
certain regions on the upper surface which may become critical during boost. 

A sketch of the MDC mated configuration is included in Figure 4. The 
NAR/GD Phase B configuration is similar. However, the orbiter nose is forward of 
the booster and the booster has a low delta wing with high canard instead of the 
high wing with low canard shown in Figure 4. As shown in Figure 5, the flow field 
produced in this mated configuration is quite complex, with interfering flow and 
unpredictable heating patterns, even at zero angle of attack. 

Examples of the flow field during reentry for the unmated configurations are 
shown in Figure 6. Because of the high booster angle of attack during 
reentry, the flow is similar to a cylinder in crossflow on the lower surface and 
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the flow over the major portion of the upper surface is separated. Locally high 
heating at the canard and below the wing still occurs. The flow is laminar during 
early reentry with transition to turbulent flow occurring slightly before peak 
heating is encountered. Thus, the heating (to lower surface areas, at least) is 
turbulent in boost, and transitional and turbulent during reentry. 

The orbiter angle of attack is lower than that of the booster for the high 
cross range 2040 km (1100 nmi) mission. Portions of the upper surface area are 
in a region of separated flow, but (depending on the geometric details) reattach- 
ment and vortex impingement may be present. Critical heating occurs on the lower 
surface and the flow is likely to be turbulent during peak heating. 

Because of the distinct differences in the heating phenomena for booster and 
or k-^ er (which differ, as well, on lower surfaces area, shielded areas* and regions 
of shock/f lowf ield interference) and the lack of complete flowfield solutions for 
such configurations, heating predictions must rely on data correlations. These 
data correlations will probably take a different form for the various configurations 
and locations and will have different levels of uncertainty. 
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III DATA COLLECTION AND EVALUATION 

The basic approach selected for this study consisted of a "hard" look at the 
ground test data generated by the Phase A and Phase B contractors. To assist in 
this evaluation , data available through SADSAC and contractor reports were cataloged 
and screened, reduced where necessary, and analyzed and correlated. 

As part of this task the tests were cataloged and summarized by configuration, 
facility, and test condition. The number of tests run in each tunnel for each 
configuration are listed in Figure 7. This chart shows that the bulk of the data 
were obtained in the Langley Variable Density Tunnel and that only a few configurations 
were tested in more than 3 facilities. It should be noted that all data were not 
available for analysis and correlation in this study, since a number of tests were 
run too late to be incorporated. In Figures 8 and 9 the Mach/Reynolds number 
simulations provided for booster and orbiter are shown. Only a few booster runs 
were made below Mach 7, a probable staging velocity for the heat sink booster. It 
should be noted, however, that higher staging velocities were under strong consider- 
ation in Phase B. The Mach simulations were even poorer for the orbiter, with no 
data taken above Mach 16 where peak heating is likely to occur. Although Figure 
9 was developed for the MDC orbiter, similar conclusions are inferred for the 
configurations of other contractors. 

Data analysis was made more burdensome by the form of the data provided by 
SADSAC, particularly when the paint technique was used to define the heat transfer. 

The heat transfer coefficients were derived for arbitrarily assumed recovery 
temperatures and nondimensionalized to reference sphere heating. The relationship 
between the q and h thus varied for the various test conditions and when compared 
with theory. The reported h had to be adjusted to the corresponding theoretical 
recovery temperature or rather large discrepancies would result. Figure 10 shows 
the relationship between q and h as a function of wall temperature and assumed 
recovery temperature ratio. 

The SADSAC isotherm data (Figure 11) did not include grids to help define the 
location of the iso therms , and a time consuming scaling of the figure was required. 
Although a corresponding grid figure was provided, stretching of the paper during 
reproduction did not allow direct overlay. 


1 1 


MCDO/V/VC(.L OOL/OL4S 45rirOA/4i;r/CS COM ¥> A MY - 



FIGURE 7 

SUMMARY OF HEAT TRANSFER TESTS ON SHUTTLE CONFIGURATIONS 


30 JUNE 1972 


MDC E0639 


P FACILITY 1 

I ARC 1 

40 

co 


CO 00 
CM rH 



rH CO 
CO CO 




AEDC 

B 

n 


1 

— 

17 




_16j 

DQ 

14 

17 

1 

40 

23 



r-* 

rH 

I CAL 

HST 

23 


1 

26 





b 

<N 

20 

20 

■ 

00 






UPWT 

CM 

CO 

11 

1 


16 




u 

oS 

(0 

*4 

CFHT HRNT 


Nf 

CO 

1 






■ 

sO UO 
CM CM 

<o 

rH r-» 
CM *cf 

21 

21 

3 

-4“ CO 
CN Mf 

29 



| VDT 

* 

CO CN 00 
CM CN 

r** Os r*> no so 

rH so r-t CO 

1 

CM CM sC 40 
CM CO CO 

tn so co m 
^ CO M CM 


OS CM CO 

CM Os 


i-4 

W 

Q 

2 

c/5 

Vi 

0 

CN 

1 

vO 

rH CN CO 40 
(Q PQ 03 CN 

0 

1 

CN 

1 CO 

CP PQ > w 

« n m tn s i oo 
CO os H H OS os H 
ill lilt 
pQ PQ OQ pQ CQ PQ OQ 

u 

8 

LaRC Clipped 
Delta 

MDAC 01 
02 
03 

050 /B 
Int Tank 
Orb 

NAR 130C 
134B 
161B 

VC70-0110D 

GRUM R0S-NB2 
H-3T 

LMSC LS-200-5 

With Ext. 

Cone/Cyl 

Tank 

LaRC 9.90 j 

VEHICLE 

BOOSTER 

ORB ITER 






12 


MCOOA/^f'LL t>OUGLAm 4SriVOM4(yriCS COMP4/VK 


NUMBER OF TEST RUNS 













































FREE STREAM REYNOLDS NUMBER BASED ON MODEL LENGTH 


30 JUNE 1972 


FIGURE 8 


MDC EOS39 



/** CDOMNBLL DOI/GL4S 4STffOW4(/riCS COMP4WV - F4ST 



30 JUNE 1972 


MDC E0639 


FIGURE 9 

COMPARISON OF GROUND TEST AND FLIGHT CONDITIONS 
MDAC ORBITERS 
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IV DATA CORRELATIONS 

The bulk of the data analysis and correlation effort was spent evaluating 
the GD/C and MDC (Martin) booster and NAR and MDC orbiter configurations. The 
correlations were of two general types: curve fits to reported h/h^^ data and 
development of factors to relate the measurements to an appropriate theory. The 
theories evaluated included the Eckert Reference Enthalpy j pp ,and Beckwith and 
Gallagher for laminar flowjand Spalding-Chi, py ,and Beckwith and Gallagher for 
turbulent flow conditions. Local flow properties were based on equivalent sharp 
cone flow except for the booster at zero degrees angle of attack for which ex- 
pansion from normal shock entropy was assumed. All theoretical calculations were 
performed with the MINIVER computer program developed by MDC and described in 
Reference 1. 
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1. Booster 
A. Windward Surfaces 

To assist in defining the local flow field properties on the windward surface 
of the GD/C boos ter , shock wave angles were measured from Schlieren photographs. 

Figure 12 shows the windward 9hock wave angle with respect to the free stream 
at the aft end of the GD/C cylindrical booster as a function of angle of attack. 

The shock angle is compared with that obtained from Delta Wing Test (Reference 2) 
and with analytical predictions for the delta, cone, and wedge. At an angle of 
attack of tt/3 radians (60 degrees) the figure show the difference between the 
shock wave and the GD/C booster surface angle to be less than tt/ 180 radians 
(1 degree), but is approximately 2 tt/ 45 radians (8 degrees) for the AFFDL delta 
wing. Although small, these differences have a significant effect on the boundary 
layer edge temperature and velocity. 

Heat transfer data obtained on the forward lower surface centerline of the 
GD/C booster were correlated with theoretical heating rates given by the 
Eckert strip, py strip, and Beckwith and Gallagher swept cylinder theories. For 
local flow deflection angles less than tt/3 radians (60 degrees) the Eckert and py 
theories were computed using the conical shock flow field assumption contained in the 
MINIVER program, and for local angles equal to, or greater than,ir/3 radians (60 degrees) 
shock wave angles were specified. Cross flow corrections, h /h . were 

computed and least square curve fits obtained. The standard deviations of 
h DATA /h THEORY about tbe curve fits were computed for each theory. The curve 
fits, their equations, and standard deviations for an angle of attack equal to 
w/9 radians (20 degrees) are shown in Figures 13 through 15. A comparison of these 
results shows that the better fits are obtained with the Eckert, and Beckwith and 
Gallagher theories. These curves show that none of the theories is able to 
model the heating, since they do not account for the three dimensional 
character of the flow. A comparison of the standard deviations throughout the 
entire angle of attack range tested, tt/9 to tt/3 radians (20 to 60 degrees), showed 
that the Beckwith and Gallagher theory yields the lowest values. However, when the 
standard deviations are compared on a percentage basis, there is little difference 
in the quality of the data fits for the three theories. 
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FIGURE 12 

SHOCK WAVE ANGLE COMPARISON 



a -ANGLE OF ATTACK ~ RADIANS (DEGREES) 
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The AEDC Tunnel B data on the GD/C B-15B-2 booster were compared directly with 
theory. A comparison of the zero angle of attack data with laminar Eckert* 
pu strip, and Beckwith and Gallagher swept cylinder results is shown in Figure 16. 
Although there is good agreement with all three theories on the extreme forward 
portion of the vehicle, little agreement exists aft of X/L = 0.4. Therefore, the 
design curve based on direct correlation in terms of h/h^ ef was used for the 


nominal case. 

A comparison of transitional and turbulent data with turbulent Spalding-Chi , 

PV ,and Beckwith and Gallagher swept cylinder results, as well as with the laminar 
theories, was also made for the data taken at high angles of attack. Typical results 

for an angle of attack equal to tt/3 radians (60 degrees) are shown in Figure 17. The 
Beckwith and Gallagher results show good agreement with the turbulent data, but the 
three dimensional character of the flow prevented the laminar theories from match- 
ing the test results. Therefore, the laminar design curves shown in this figure 
are based on direct correlation of the data. The ratio of data to theory, 

h /h . for the turbulent Beckwith and Gallagher results at tt/3 radian 

DATA THEORY* 

(60 degree) angle of attack is shown in Figure 18. This correlation was used for 

the nominal turbulent, tt/3 radians (60 degree) angle of attack case. 

Peripheral heating data for both the straight and delta wing booster configur 

ations at X/L =0.25 and angle of attack equal to tt/6 radians (30 degrees) are shown in 

Figure 19. These data were obtained from the NASA Langley Continuous Flow Hypersonic 

and Variable Density Tunnels and the NASA Ames 3.5 foot tunnel. The data for these 

various configurations have been combined because the geometric similarity in the 

forward portion would be expected to produce equivalent peripheral distributions. 

This figure also shows a least squares curve fit to the data in the form 

In (h/'h ) = f(0, 0 2 , 0 3 ) with the angle 0 measured from the lower surface center- 

ref 

line. In general, the Mach 7.95 paint data are higher than either the Mach 10 paint 
and gage data or the Mach 7.4 gage data. Since no reasonable explanation of these 
differences could be determined, all of the data were used in the correlation. 

General trends of heat transfer data on the wing lower surface of the GD/C 
B-15B-2 booster were investigated by CalComp plotting from a data tape obtained 
from the AEDC Tunnel B facility. Machine plots of the wing lower surface heat 
transfer coefficient, as a function of percent span, were made for eight chordwise 
locations over the angle-of-at tack range tested. From these plots it was determined 
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FIGURE 16 

GD/C (B-15B-2) BOOSTER LOWER SURFACE CENTERLINE HEATING DISTRIBUTION 



25 


MC0o/v«n.L actual. Am Am-rmo*AUTicm company 


KAmr 















Ref 


30 JUNE 1972 


MDC E0639 


FIGURE 17 

GD/C B-15B-2 BOOSTER LOWER SURFACE CENTERLINE HEATING DISTRIBUTION 
ANGLE OF ATTACK = — %r~ RADIANS (60°) 
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FIGURE 18 

RATIO OF LOWER SURFACE CENTERLINE HEAT TRANSFER COEFFICIENT DATA 
TO TURBULENT BECKWITH AND GALLAGHER SWEPT CYLINDER VALUE 
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FIGURE 19 

PERIPHERAL HEATING DISTRIBm'ION 

ANGLE OF ATTACK = RADIANS (30°) 
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X/L = 0.25 



28 

MCDOWWFLL DOI/OL4S ASTWOWAt/TICS COMP/XAIY - E4IT 








30 JUNE 1972 


MDC E0639 


that as the angle of attack increases to tt/3 radians (60 degrees) the heat transfer 
coefficients tend to become relatively independent of chordwise and spanwise loca- 
tion except in the vicinity of the leading edge. 

Data-theory comparisons were made in the chordwise direction at the 50 percent ^ *a 
location on the lower surface. Eckert laminar, and Spalding and Chi turbulent strip 
theory, calculations were made using both oblique and conical shock flow field 
assumptions. These calculations were carried out over the total surface distance 
at this spanwise location. Also, in conjunction with the Eckert, and Spalding and 
Chi calculations, a Prandtl-Meyer expansion was carried out from the 50 percent chord 
location to the trailing edge of the wing. There was little difference in the 
results regardless of method, but the oblique shock calculation showed slightly 
better agreement with data at the lower angles of attack. On the basis of the 
theoretical calculations it was concluded that the flow was laminar for angles of 
attack of 0 and tt/18 radians (0 and 10 degrees), turbulent at tt/9 radians (40 degrees) 
when grit was applied to the lower surface and laminar without the grit. The 
5 tt/ 18 radian (50 degree) data appeared to be transitional and the tt/ 3 radian (60 degree) 
data were turbulent (with or without grit on the lower surface) . The data and 
their respective curve fits used for the nominal case are shown at the 50 percent 
span location for zero and tt/3 radian (60 degree) angles of attack in Figures 20 and 21. 
B. Leeward Surfaces 

Upper surface centerline data for the GD/C B-15B—2 booster at zero angle of 
attack, obtained in AEDC Tunnel B, are shown in Figure 22. The curve fit shown in 
the figure was used in conjunction with mated interference results as the nominal 
case . 

Upper surface centerline data on the GD/C B-15B-2 at tt/3 radian (60 degree) 
angle of attack are shown in Figure 23. A portion of these data were obtained 
during a test with grit on the lower surface to induce turbulent flow. As indicated 
by the data, the grit did not affect the upper surface centerline flow. These 
data vary around the Eckert laminar strip theory computed for zero angle of attack. 

The systematic variations in the data indicate that this is not data scatter and 
is probably due to vortex impingement on the upper surface centerline. Because 
the Eckert laminar theory went through the mean of the data, it was used in the 
estimates of the thermal protection system weight. 

Wing upper surface heat transfer data at the 50-percent span location is shown 
plotted against angle of attack in Figure 24. At the lower angles of attack the 
data decrease moving aft along the chord, while at the higher angles of attack 
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FIGURE 20 

GD/C B-15B-2 BOOSTER WING LOWER SURFACE HEATING DISTRIBUTION 
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FIGURE 21 

GD/C B-15B-2 BOOSTER WING LOWER SURFACE HEATING DISTRIBUTION 

ANGLE OF ATTACK: a - JL. RADIANS (60°) 

Moo - 8 .0 , 3 

Reooftn - 12.3 x 10 (3.75 x 10VFT) 

SEMI-SPAN LOCATION: Y/S - 0.5 
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FIGURE 22 

GU/C B-15B-2 UPPER SURFACE CENTERLINE HEATING DISTRIBUTION 
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FIGURE 23 

GD/C B-15B-2 UPPER SURFACE CENTERLINE HEATING DISTRIBUTION 
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FIGURE 24 

GD/C B-15B-2 BOOSTER WING UPPER SURFACE HEATING DISTRIBUTION 

M = 8.0 

oo 
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when the upper surface becomes shielded there is very little difference at the 
three chordwise locations. No systematic correlation of these data with the 
Eckert strip theory was found, and the data fairings shown on the figure were 
selected as nominal. 

C. Interference Regions 

The portions of the booster likely to be affected by interference are in the 
vicinity of the canard, the aft side body above the wing, and the upper surfaces 
during mated boost. 

At zero angle of attack the mated interference overpowers any canard effects. 
Figure 25 shows the peripheral distribution of the heat transfer coefficient in the 
canard region at tt/3 radian (60 degree) angle of attack. Some of the thermocouples 
were covered by the canard, and the data presented in the SADSAC report from those 
locations are not valid. These have been identified in the figures. The dashed 
line shown on the figure was generated by using a mean peripheral heat transfer 
distribution curve (h0/hg) from Reference 3 for a clean cylinder at Mach 8.0. This 
was then applied to the centerline value of h/h^^. This figure shows that the 
canard has no discernible effect on the heating distribution. The probable reason 
is that the canard is at zero angle of attack relative to the tunnel flow and the 
region of the body likely to be influenced by the canard is in a region of 
separated flow at this angle of attack. 

The ratios of the local heat transfer coefficients on the side body above the 
wing (measured on the B-9J and B-15B-2 boosters) to those on the lower surface 
centerline were determined at several axial locations. The peak local values are 
presented in Figure 26. The trend of the data shows a decrease in peak value of this 
ratio with increasing angle of attack. This is due to increases in the lower sur- 
face heating with angle of attack, rather than to decreasing values of the body heating. 
For zero angle of attack the multipliers vary between 2.6 and 9 times the lower 
centerline value. 

Interference heating on the mated GD/C booster was evaluated using the data of 
Reference 6. The configuration was comprised of the GD/C B-15B-2 booster and the 
NAR 161B delta wing orbiter. The ratios of the local heat transfer coefficients 
on the mated booster upper body to those on the unmated booster were determined 
as a function of axial location for four peripheral angles. Typical results are 
shown in Figures 27 and 28. As may be seen in the figures, heat transfer 
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GD/C B-15B-2 BOOSTER PERIPHERAL HEATING 
CANARD REGION 
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FIGURE 26 
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coefficients on the mated booster in some locations are as much as 4 times those 
on the booster alone at 0 = it/ 2 radians (90 degree) increasing to approximately 20 
on the top of the vehicle. The data show no systematic trend with Reynolds number, 

A correlation of mated interference effects as a function of free stream 
Mach number, obtained from Reference 8, is shown in Figure 29, Superimposed 
on the data of the reference are the two peak values on the upper centerline 
of the B-15B-2 booster (Figure 28), These data show very good agreement with 
the correlation and were therefore used to determine the nominal mated inter- 
ference effects. The correlation assumes a Mach number dependence on the ratio 
of mated to unmated heating. Alternate correlations involving Reynolds number 
were also evaluated, but provided no better fit. 

Upper wing data on the mated B— 15B— 2 booster exhibited more scatter than on 
the unmated configuration, however, no net increase of the heating level was 
noted. In addition, no bow shock interactions were observed on the instrumented 
portion of the wing which extended outboard to 70 percent of the semispan. 
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FIGURE 29 


CORRELATION OF INTERFERENCE HEATING DATA 
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2 Orbiter 

a. Windward Surfaces - Lower surface centerline data correlations included 
both least-square curve fits and multiple regression analysis of the NAR and 
MDC orbiter data. 

Least-square curve fits were made to the ratio of the Reference 9 data to 
theory on the NAR orbiter lower surface centerline for the pp i and Eckert laminar 
strip theories and Beckwith and Gallagher laminar swept cylinder theories. 

Resultant curve fits at an angle of attack equal to ^ / 6 radians (30 degrees) for 
the pp and Beckwith and Gallagher theories are shown in Figures 30 and 31, 
respectively. Similar fits were made to the MDC delta wing orbiter data of 
Reference 10. Standard deviations of the data about the curve fits were obtained 
and typical results are shown for the MDC orbiter in Figure 32. As seen in 
this figure, there is little difference between the three theories, except for 
Beckwith and Gallagher at zero angle of attack where this theory is not expected 
to be valid. 

Multiple regression analysis was employed to refine the computed 

data fits. Given a set of observations of the dependent and independent 

variables involved in a particular experiment, multiple regression can be 

employed to obtain a best fit to the data by an equation of the form: 

y = a + a.x- + a a x n + a~x_ + . . . + a x 
J o 11 22 33 nn 

where y is the dependent variable and x, x^ , « • • » x n are independent vari- 

ables. A multiple regression solution gives the least-square best value of 

the coefficients a , a, , . . . , a . In addition to the final equation, inter- 
o 1 n 

mediate regression equations are obtained after each step in the MRA which give 
an indication of which variables are most important to the analysis. Also some 
of the variables included in the analysis may be rejected because they have no 
significant effect on the dependent variable. For each step in the MRA, pre- 
dictions of the dependent variable (based on the current equation) are determined 
along with statistical information regarding goodness of fit, multiple correlation 
coefficients , and various measures of the reliability of the coefficients. Of 
particular importance is the standard error of estimate for each step which 
represents the RMS error of the predictions. MRA is a powerful tool when 
little is known about the functional form of the dependent variable, but the 
method is enhanced when such knowledge is available. A more thorough 
discussion of MRA is presented in Reference 11. The functional 
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FIGURE 30 

LOWER SURFACE CENTERLINE DATA CORRELATION FOR NAR DELTA WING ORBITER 

Angle of Attack: “ = w /6 RADIANS (30°) 

Theory: PH Laminar Strip Theory 
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co 

Reference: 9 



MCDOWM1.L DOUOLJk9 4STJV0M40rfCS o 




DATA THEORY 


30 JUNE 1972 


MDC E0639 


FIGURE 31 

LOWER SURFACE CENTERLINE DATA CORRELATION FOR NAR DELTA WING ORBITER 
Angle of Attack: a = tt/6 RADIANS (30°) 

Theory: Beckwith & Gallagher Laminar Strip Theory 

Flagged Symbols: M„ 10.0 

Unflagged Symbols: =8.6 

Reference: 9 
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FIGURE 32 


COMPARISON OF STANDARD DEVIATION FOR DIVERGENCE PARAMETER 
MDAC DELTA WING ORBITER 
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form assumed for the NAR orbiter data was h_ /h_, = f [X/L, (X/L) , M a > 

Data Theory 

ReQQ /L ] • Results obtained are shown in Figure 33. The standard deviation is 
reduced with the use of this method but not appreciably. 

Correlation of the MDC lower surface centerline data included least-square 
curve fits to the ratio of data to theory. The curve fits, shown in Figure 34 
were obtained using the Eckert theory for laminar flow and the Spalding and Chi 
theory for turbulent flow. Crossflow corrections as described in Reference 12 
were applied to the theories. The curves shown are third order fits to the local 
flow deflection angle, 8. Multiple regression analysis was also applied to the 
MDC orbiter data obtained at AEDC Tunnel B and Cornell. Three cases considered 
are summarized in Figure 35. The stepwise development of the MRA equation is 
shown for each case. The first case employed flow deflection angle, velocity 
and unit Reynolds number as independent variables. The second case dropped the 
variable 6 and the third case employed powers of <5 up to third order and, hence, 
gives the same results as the least-square fits. As shown in the figure, Case 1 
yields a lower standard error and a somewhat better fit than Case 3, as indicated 
by the correlation coefficient. As with the NAR data, the differences are small 
and it was concluded that the use of least-square curve fits was the more expedient 
method of obtaining results within the accuracies required. Therefore, the curves 
of Figure 34 were used as nominal for TPS design. 
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FIGURE 33 

COMPARISON OF LEAST SQUARE CURVE FIT WITH MULTIPLE REGRESSION ANALYSIS 

NAR DELTA WING ORBITER 
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FIGURE 34 

MDAC 050/B ORBITER LOWER SURFACE CENTERLINE DATA CORRELATION 
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FIGURE 35 

STEPWISE MULTIPLE REGRESSION ANALYSIS 
OF MDAC/DWO LOWER SURFACE CENTERLINE DATA 
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Lower surface spanwise heat transfer distributions on the MDC orbiter were 
determined utilizing Shuttle Phase B results for both the MDC and NAR orbiters. 
Nominal spanwise distributions normalized by the centerline laminar heat transfer 
coefficient, obtained from Reference 13* are shown versus angle of attack for 
three spanwise locations at X/L = 0.5 in Figure 36. Similar results are shown 
in Figure 37 for an angle of attack equal to 0.436 radians (25 degrees) at X/L ■ 0.75. 
A comparison of the data scatter on the MDC orbiter with that on the NAR orbiter, 
obtained utilizing the data of Reference 9, shows good agreement. The nominal 
spanwise distributions with turbulent centerline flow, also obtained from 
Reference 13, are shown against angle of attack for three spanwise locations 
in Figure 38. 

Upper Surfaces - Upper surface centerline data for the MDC orbiter, 
obtained in AEDC Tunnels B and F, are shown in Figures 39 and 40 for angles of 
attack equal to 77/ 6 and 0.698 radians (30 and 40 degrees). As shown in the figures, 
separate curve fits were used forward and aft of the canopy because of the large 
difference in heating levels. Although there is little difference in the results 
at the two angles of attack, a design value for the nominal case was interpolated 
for the 0.593 radian (34 degree) angle of attack. 

The upper wing data for the MDC orbiter obtained at Cornell were compared 
with Eckert laminar strip theory. Such a comparison was made since one traditional 
method of predicting heating rates on leeward surfaces in separated flow has been 
the use of a percentage (50 to 55 percent) of flat plate heating rates. The data 
exhibited a scatter of one order of magnitude, as shown in Figure 41, and the 
highest data points were less than 20 percent of the theoretical results. Also, no 
systematic variations with Mach number or Reynolds number could be found. It was, 
therefore, concluded that the linear fit shown on the figure would be the most 
reasonable for the nominal case. 

In order to obtain a better understanding of the leeward side flow character- 
istics on the MDC orbiter, oil flow data from tests in the Ames 3.5 ft hypersonic 
tunnel were studied. The data show that flow separation occurs as the flow 
rounds the leading edge chine for angles of attack from 0.35 to 0.87 radians 
(20 and 50 degrees). Location °f the separation varied with angle of attack as shown 
in Figure 42. These results confirm, the low heating levels obtained on the 
upper wing. The data also indicated a weak vortex flow which streams along the 
upper body centerline and impinges on the canopy. The relative weakness of this 
impingement phenomenon was verified by the Cornell tests on this configuration. 

50 


MCDOMA/fl.1. DOUGLAS 49TJ70A/4L/T/CS COMP4MK - 8 T 



30 JUNE 1972 


FIGURE 36 


MDC E0639 



0 IT/ 18 7T/9 7T/6 27T/9 57T/18 7T/ 3 

(0°) (10°) (20°) (30°) (40°) (50°) (60°; 

ANGLE OF ATTACK - RADIANS (DEGREES) 


51 


Mcoow^ffu OOOG14S 4sr»ow>u;r/c5 co/wp>i/vv - 






30 JUNE 1972 


MDC E0639 


FIGURE 37 

LOWER SURFACE SPANWISE HEATING DISTRIBUTION 
FOR MDAC DELTA WING ORBITER (050/B) 
LAMINAR CENTERLINE FLOW, X/L = 0.75 
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FIGURE 38 

MDAC DELTA WING ORBITER (050/B) LOWER SURFACE SPANWISE HEATING DISTRIBUTION 
TURBULENT CENTERLINE FLOW, X/L = 0.5 
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3. Boundary Layer Transition - Several recent studies of the problem 
of transition during hypersonic entry have focused on parameters derived from 
sophisticated flow field descriptions. The availability of electronic computers 
has allowed application of this approach to axisymmetric flow fields (Reference 
14). For more complex flow fields, complete solutions are not possible, and 
restrictive assumptions are necessary (Reference 15). The correlation of data with 
relatively simple parameters which may be regarded as indices of the complex 
phenomena was applied for this study of Shuttle transition design criteria. The 
response of a flow field to variations of a simple parameter may, however, be much 
more pronounced for one configuration or flow environment than f° r another. As an 
extreme example, angle of attack can drastically alter the flow field about a cone, 
but is an irrelevant parameter for the flow about a sphere. Therefore, correlations 
based on simplified parameters are likely to be applicable only to the class of 
configurations and attitudes represented by the data evaluated. 

The Phases A and B boundary layer transition data have been recently 
generated, and no comprehensive survey and correlation study of those data were 
available prior to this program. This survey, analysis , and correlation identified 
a group of 29 data points which were felt to be reliable. These data, which were 
gathered in continuous flow, shock tunnel , and hot shot facilities using (in most 
cases) delta wing orbiter models, are summarized in Figures 43 and 44. Also 
included are most of the basic correlating parameters used during this study. Local 
flow conditions assumed a conical flow field. Plots of these data indicated 
regular variations with several commonly applied transition parameters. Correlations 
in terms of local and momentum thickness Reynolds numbers are shown in Figures 45 
and 46. Comparison of the figures show that the booster transition data correlate 
better with the orbiter data if the Reynolds number based on boundary layer 
momentum thickness at transition onset is used. Figure 47 shows the same data 
correlated in terms of the MDC Phase B transition criterion. Included in this 
figure are the criterion and a least-square fit to the data. It should be 
noted that the Phase B criterion was a fit to a large number of ground and flight 
test data on delta configurations, whereas the data in the figure were all 
obtained on Shuttle shapes. For purposes of this study the least square fit was 
used as the best fit correlation. 

Evaluation of the data in each of the above correlations shows that the mean 
of the transition data taken in the continuous flow AEDC Tunnel B is higher than 
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FIGURE 45 
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MDAC PHASE B TRANSITION CORRELATION PARAMETER 
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the mean of the data taken either in the Cornell Hypersonic Shock Tunnel 
or AEDC Tunnel F. These differences were found to be less pronounced if a wall 
cooling parameter was introduced. In particular, the use of a parameter T^/T^ 
was found to reduce these tunnel differences. Alternate correlation f rms were 
evaluated as a part of I RAD at MDAC-E, and are included herein for information. 

The same data points were processed by multiple regression analysis computer 
program assuming various combinations of parameters and forms of correlating 
equations. Those multiple regression fits which provided a significant degree of 
correlation are presented in Figures 48 through 50 along with their one-standard- 
error-of-es timate deviations. It should be noted that momentum thickness Reynolds 
number is roughly proportional to the square root of transition length (this is 
not a direct proportionality due to application of streamline divergence factors 
calculated by the method of Reference 16) . Thus, momentum thickness Reynolds number 
is less sensitive to transition length than either length Reynolds number or 
fractional body length (X/L) . As a result, equations based on momentum thickness 
Reynolds number exhibit generally lower standard errors of estimate than equations 
with the other dependent variables. For a given dependent variable, however, the 
standard errors of estimate gave a quantitative basis for ranking the equations. 

Equations including Reynolds number, Mach number, flow deflection angle, and 
either the ratio of wall temperature to static temperature or the ratio of wall 
temperature to total temperature were found to yield low standard errors of estimate 
and probability plots consistent with the assumption of normal distributions. A 
minor improvement in the standard error of estimate was achieved by inclusion of 
unit Reynolds number, but the indicated dependency is slight. Equations including 
a cubic expression in flow deflection angle were also effective, but such com- 
plexity often yields misleading results when applied to a small data sample. 

The significance of the correlations to vehicle design must be evaluated 
in terms of flight parameters. Thus, studies discussed in this section can only 
be used for the screening of correlations. Altitude uncertainties associated with 
the correlations are a better measure of the degree of correlation provided. 

Results of this approach are included in Section V of this report. 
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FIGURE 48 

UNCERTAINTY IN TRANSITON ONSET BASED ON LOCAL REYNOLDS NUMBER CORRELATIONS 



65 


mcoowfll oot/GL4s 4sr«vo^4i;ncs 


COMP4MK - f4tr 



30 JUNE 1972 


FIGURE 49 

UNCERTAINTY IN TRANSITION ONSET BASED ON 
MOMENTUM THICKNESS REYNOLDS NUMBER CORRELATIONS 
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FIGURE 50 


UNCERTAINTY IN TRANSITION ONSET LOCATION 
BASED ON COMBINED REYNOLDS NUMBER AND ANGLE OF ATTACK CORRELATIONS 
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V. DEFINITION OF AERODYNAMIC HEATING UNCERTAINTIES 

As shown in the correlation curves of Section IV, the analyst is faced with 
aerodynamic heating data which scatter significantly about theory or the selected 
correlations. Therefore, he must decide on factors to apply to his best fits which 
assure a safe design for flight. In this study an approach was developed and applied 
to quantify the uncertainties for various confidence levels using a statistical 
approach. This approach is modeled after procedures utilized in the Dynasoar program 
and successfully applied to the design of the hypersonic BGRV flight vehicle. 

In the BGRV program ground test data were correlated in terms of significant 
aerothermodynami c parameters. Best fits to the data were then made and the standard 
deviations of the data about the best fits were determined. Factors were developed 
to apply to the nominal data fit to estimate the off-nominal effects of aero- 
dynamic heating on temperature and TPS weight. The increments (in temperature and 
weight) were then statistically combined with other independent factors which affect 
the temperature and TPS weight using the root sum square method. These included, 
for example, trajectory dispersions, off-nominal atmospheric properties, and varia- 
tions in material properties affecting TPS performance (surface emittance, skin 
thickness, skin and insulation properties). 

Post flight evaluation of the BGRV data showed that the data scatter about the 
ground test developed correlations were essentially bounded in terms of peak heating 
by the band representing two standard deviations of the ground test data about its 
mean. These results provide support for the use of ground test data and data 
scatter to establish uncertainties in flight even though some of the factors which 
produce the ground test scatter do not occur in flight. It should be recognized, 
however, that the BGRV was a relatively simple (axi symmetric) shape with a high 
quality ground test data base which made the uncertainty low, compared with the 
currently available Shuttle data base. Thus,direct application of the results of 
this study are expected to be conservative. However, as the Shuttle Phase C/D 
data base is developed, it is anticipated that better estimates of the effect of 
uncertainties in the heating prediction methods can be obtained with this 
statistical approach. 

1. The Statistical Approach - Since it is suggested that a statistical 
approach be used to develop heating factors for flight, it is necessary to select 
and apply correctly, statistical methods which model the problem. 
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The specific problem is to apply a factor to heating analysis methods that will 
assure, with some level of confidence, that the predicted heating will not be exceeded 
in flight for known free-stream conditions. The prediction methods which are 
currently considered state of the art are not exact for Shuttle-like configurations 
and it is necessary to rely on correlations developed from comparisons of ground test 
measurements with these imperfect theories. With this approach a factor, or 
factors, is determined to apply to the theory in order to scale to flight. Since 
we are dealing with data there is scatter. This may be due to tunnel measurement 
errors or to an inexact match of the theory to the physical problem. If the latter 
case is assumed, the ground test scatter and deviation about the data mean can be 
used to define a scatter or uncertainty band for flight. To assist in 
selecting the appropriate statistics to apply, a few definitions and statistical 
concepts are described below. 

a. The Normal Distribution - The normal (Gaussian) probability distribution 

2 

is completely determined by its mean, /x,and its variance,^ . The statement that 

2 

a variable, x, is normally distributed with mean p and variance cr can be written 

2 

x : N(/x,a ) 

Tables of the normal distribution are available which provide tabulated values 
of K, the probability, P, that a randomly occurring value of the variable, x, will 
be within an interval of length, Ka, on either side of the mean , /x . That is, the 
probability that 


/U - Ka < x < + Ka . 

b. Normal Distribution Sample Statistics - A sample of size n from a normal 
distribution is a collection ( x 2 ,,X n^ norma lly distributed variables. The 
sample mean is given by 


~ l n 

X = - I x . 

n i=l 1 


and the sample variance, s , is given by 


( 1 ) 


9 1 n — ? 

s = T I (x-x ) - (2) 

i=l 

When drawing inferences about the normally distributed parent population (or 
universe) from which the sample was taken, we use x and an unbiased sample 
estimate of the population variance , , as defined by 


*2 = 


n 

I 

i=l 


- x ) 2 


( 3 ) 
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It is important to recognize that x. is not u and o 2 is not o 2 since the sample size 
is limited and does not include the total population of data. However, these 
parameters and the assumption that we are dealing with normally distributed data 
can be used to predict the characteristics of other data samples from the total 
data population. 

An important property of the normal distribution is that a sample mean, 7, 
for a random sample of size n is normally distributed about the true mean p with a 
variance a /n. Another important property is that the quantity ns 2 /a 2 has an x 2 
(Chi-square) distribution with (n-1) degrees of freedom, where from equation (2) 


2 n 

n s = S (x. 
i-1 



(4) 


The Chi square variable has (n-1) degrees of freedom rather than n because 

only (n-1) of the x ± are "free", the other one being determined from the fact that 
x is known and 


n x 


n 


£ 

i-1 


x . 


l 


(5) 


c * -Confidence Limits for Normal Population Parameters — A major use of 
statistics is to use sample data to draw inferences about the parent population. 
Limits can be specified which include the population mean and variance within a 
stated level of confidence if the data are normally distributed. In this study, 
the sample data are derived from ground tests and are assumed to be a part of the 
total population which includes ground and flight data. 

Confidence Limits on the Mean - It is possible to place confidence limits 
on the mean of the total population based on the sample measurement, e.g., to 
determine limits on the mean in flight based on the ground sample. To do this, 
it is necessary to recognize that the variable 

t (6 , 

has a Student's t distribution (reference 14) with (n-1) degrees of freedom. 

2 a 9 

Because s and o+ 9 are related by the expression 

o ' 2 =-EL_ s 2 ( 7 ) 

n-1 


equation (6) can also be written 
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. _ (x - y) 
t 7C 




where 8 is defined, for this purpose, by the equation o = 0 ^ 

Using equation (8) , we can define an interval which includes the mean based 
on the sample data. For example, if t Q represents the value of Student's t 
corresponding to a probability of 0.05 that |tj > q^ 5 there is a 95-percent 


chance that 


( x - /*) fo L L c (9) 

a P C 0 , 05 

a 1 

where t 0.05 is obtained from tabulated results for Student’s t for (n-1) degrees 

of freedom. Thus, a 95-percent confidence interval for the population mean is 

defined by the expression A 

g_ u ^ - a (10) 

x - *0.05 vPT < x+ *0.05 -fzr- 

If some probability other than 0.05 is desired we replace t^ by the corresponding 
t in the above expression Equation (10) can also be used to define a confidence 
interval for the ratio of the population mean to the sample mean by rewriting it in 
the form 

jl = i + k _<L an 

x “ x n nr 

Confidence Limits on the Variance - Confidence limits on the variance 

9 ~~~ “ 2 2 2 ... . 

( *“) can be derived from the fact that ns / c r has a X distribution with 

(n-1) degrees of freedom. If the confidence level is to be 0.95, that is, if an 

interval is desired such that in repeated experiments 95-percent of the intervals 

2 

derived in the same manner would contain the true population variance, o , then 
the interval for a particular sample of size n would be derived from the 


re 1 at ion , 


2 < n s' 


where X 2 and X 4 - _ are derived from Chi-square tables for (n -1) degrees of freedom 
12 2 l 2 2 

no that the probability is 0.975 that X > X^ and 0.025 that X _£ X ^ 

One-Sided Confidence Limits - The above discussion assumed that both 
upper and lower confidence limits for the mean and variance were of interest. if, 
as is the. case for heating uncertainties, we are really interested that a heating 
rate or value not be exceeded at some confidence level, we are dealing in one-sided 
confidence limits, i.e., we are interested in the precent of samples below some limit 
rather than between two equally spaced limits. The symmetric nature of the normal 
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distribution curve allows calculation of one-sided confidence limits on the mean 
by evaluating Student’s t for a probability twice the two-sided value. For 

example, if a one-sided upper 95-percent confidence interval for the mean is desired 
the value of Student's t for (n-1) degrees of freedom would be used such that the 
probability is 0.10 that Jt| > t Q 1Q . Then the probability is 0.95 that 


or that 



< t o.io 


(13) 


x + t 


0 . 10 'fn 


(14) 


The value of t Q 1Q such that Jt| ^ t Q 1Q corresponds to a probability that 5-percent 
of the Student's t distribution is In each tail of the normal distribution curve. It 

is the correct value of Student's t for a one-side (upper or lower) 95-percent 
confidence interval for \i . 


d ‘ Tole ranee Limits With C o nfidence Coefficients - The previous paragraph 
discussed in general the relations between sample data and a total population 
and ways of placing confidence limits on the mean and variance of the total 
population of which the data sample is a part. It does not, however, tell where 
the next data point is likely to fall, which is the real problem to be addressed 
in developing heating uncertainties. That is, it is desired to apply a factor 
to the average heating data that will assure that some percentage of the data 
population will not exceed the resulting off-nominal heating rate. The 
appropriate statistics to use in developing such factors must consider tolerance 
limits with confidence. If the normal distribution parameters P and are known 
for the total population we can find the appropriate values of K to establish the 
desired tolerance limits on x so that 


yu. - Kcr < x < y +K a ( 15 ) 

* 2 

In practice /i andCT are not known and tolerance limits must be estimated 
from the sample mean, x, and the unbiased sample estimate of the 
population variance, o 2 . This problem has been evaluated by 

A. H. Bowker and tables of tolerance limits with confidence coefficients have been 
tabulated and published in References 17 and 18. Bowker' s table gives values of K 
such that in a large series of samples from a normally distributed total population 

72 


MCOOWMLL OOUG143 4ST4V0A/40n0S 


- C/lflr 



30 JUNE 1972 


MDC E0639 


a certain proportion of the intervals x + K will include at least a fixed per“ 
centage of the total data. In simpler terms the approach defines the percentage 
of data which falls within some tolerance band for a given confidence level. As dis- 
cussed in Paragraph 1C, one-sided confidence limits are more applicable to the pro- 
blem of heating uncertainties. 
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2 ' °£ yelopment o f the Standard Deviation - The standard deviation of the 
data about best fits to heating data was obtained for several locations on the 
booster and orbiter configurations to assist in selection of the theory which 
best matched the data, and to define uncertainty factors to apply for the flight 
trajectory. The configurations were divided into areas for which a given predic- 
tion method (and corresponding uncertainty) was expected to apply. These 
included (1) lower surface centerline, (2) of f- centerline windward regions, 

(3) shielded regions, (4) upper surface centerline, and (5) regions of shock 
impingement and interaction. In addition, the location of boundary layer trans- 
ition on the lower surface centerline was correlated and the standard deviation 
of the data around the best fits determined. For the lower surfaces, turbulent 
and laminar data were separately correlated and the corresponding standard 
deviations determined. 

The standard deviation about least square fits of the data (or ratio of data 
to theory) was obtained graphically. The graphical technique, an example of 
which is shown in Figure 51, provides a visual description of the data scatter 
distributions and a consistent technique for deriving the standard deviation of 
the data set(s). This example curve, which is based on the orbiter lower surface 
centerline correlations of Figure 34 of Section IV, contains the ratio of data 
to theory as a function of the percent of the data which lies below each value. 
Fits to the laminar and turbulent data are superimposed in the figure. It can be 
seen that the uncertainty in the laminar data is slightly greater for the 
laminar than for the turbulent data. They are, however, sufficiently close that a 
single value could be used for each set. This value of the factor corresponding 
to one standard deviation is approximately 18 percent. The laminar theories used 
to develop this curve were Eckert’s Reference Enthalpy Method for laminar and 
the Spaldmg-Chi correlations for turbulent flow. The data were obtained from 
tests of the MDC 050 B orbiter configuration in the AEDC Tunnels B and F and the 
Cornell Aeronautical Laboratory Hypersonic Shock Tunnel. The Mach number range 
was from 7.6 to 12.5 and Reynolds numbers based on free stream conditions and 
model length varied between 0.25 and 40 million. 

The slightly higher value for laminar flow had not been anticipated, since 
laminar theories are further advanced than those for turbulent flow (at least 
for stagnation point heating). Investigation of the reasons for this has 
indicated that the range of environmental parameters was wider for the laminar 
than for the turbulent data and, in addition, that the laminar heating is more 
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sensitive to three-dimensional flow effects, which are not perfectly modeled in 
theory. Further evaluation of this phenomenon will require a more complete data 
base or the development of more sophisticated flow field and heating models 
for the complex shuttle configuration. 

The B-15B-2 booster laminar heating was correlated as a function of h/h 

T6 f 

This type correlation is extremely sensitive to the configuration. Sufficient 
data were not available on the B-15B-2 configuration to establish uncertainties 
which could be applied to flight. Therefore, data correlations of similar 
configurations were used to establish the standard deviation. The data used to 
establish standard deviations for the lower surface centerline are shown in 
Figures 52 and 53. These data were obtained on three booster configurations: 
the B-15B-2, B-9J, and B-8B, all correctable in terms of h/h ref since the nose 
geometry (forward of X/L = 0.4) was essentially the same for 111 configurations 
except for the B-15B-2 which differed aft of X/L = 0.20. This can be seen by 
reviewing the data in the figure which show little geometric sensitivity forward 
of X/L =0.04. These laminar data were obtained over a Reynolds number range 
between 0.8 and 10 million based on free-stream conditions and model length. During 
the time span in which these uncertainties were developed, the booster data were 
available in a limited Mach number range, M = 7.4 to 10. It would be desirable 
to include additional low Mach number data since the Mach numbers of interest for 
heat sink boosters are below Mach 7. The curve used to define the heating factor 
corresponding to one standard deviation is included in Figure 54 for the zero 
degree angle of attack data. By comparing this curve with Figure 51 similar 
characteristics are noted but slightly lower uncertainties are indicated for 
these booster data. This might be expected since the booster nose section is 
cylindrical and thus more easily modeled with Fay-Riddell stagnation point theory. 

Fewer turbulent data were available for the booster than for the orbiter configurations 
The only configuration for which a large sample of turbulent data were available 
was the MDC booster. These data were best fit with the Beckwith and Gallagher 
theory at 60 degrees angle of attack, and uncertainty factors developed for fits 
to this theory. The range of tunnel variables encompassed Mach numbers between 
8.0 and 11.2 and free stream Reynolds numbers based on model length between 5.0 
and 15 million. The ratio of wall to total temperature ratio for these data 
varied between 0.2 and 0.6. The heating multiplier for these data was found to 
be 1.25 for one standard deviation. 
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FIGURE 52 

LOWER SURFACE CENTERLINE HEATING DISTRIBUTIONS FOR GD/C BOOSTERS 
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LOWER SURFACE CENTERLINE HEATING DISTRIBUTIONS FOR GD/C BOOSTERS 

ANGLE OF ATTACK: a = 7T/3 RADIANS (60°) 
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PROBABILITY PLOT FOR DETERMINING STANDARD DEVIATION FOR GD/C BOOSTER CONFIGURATION 
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Heating multipliers developed in a similar manner for the various locations 
on the booster and orbiter are shown in Figures 55 and 56. These tables show that 
the largest uncertainty occurs in regions above the wing at high angle of attack. 

This is due to two factors: (1) the complex phenomena related to separated and 
reattaching flows which are poorly modeled in current correlations , and (2) the low 
heating levels in tests which are near the sensitivity level of the instro.entation. 
Also, as expected, the regions where interference heating is present (upper 
surfaces of the booster at a = 0) have high multipliers corresponding to one 
standard deviation of th$ data. The data scatter is slightly greater for the 
orbiter than for the booster. This is probably due to the greater complexity 
of the flow field (during reentry). Insufficient spanwise heating data were 
available on the lower surface of the 050/B orbiter to define a heating factor. 

Thus the factors for spanwise heating were defined from tests of the NAR 
13AB delta orbiter configuration. 

In addition to heating factors applied to laminar and turbulent heating 
methods, an uncertainty exists in the state of the flow. This is best quantified 
as an uncertainty in the time of transition from laminar to turbulent flow, which 
can constitute a major contribution • to the total uncertainty since turbulent 
heating is a factor of two or more greater than the corresponding laminar values. 

The transition data were treated similarly to the heat transfer 
results. Best fits to a number of correlating approaches were made, and the standard 
deviation of the data about the best fit determined graphically. The probability 
plot developed for correlating forms Jte Q vs M g , is included in Figure 57. The 
curve is nonlinear, suggesting that the data are not normally distributed about 
the best fit and a rather large standard deviation is noted. This is not surpris- 
ing when one considers the complexity of the phenomena and the acknowledged 
sensitivity of the transition process to wind tunnel disturbances. 

The times at which transition is predicted to occur for three of the more common 
correlations (and uncertainties corresponding to one, two and three standard deviations 
about the best fit) are summarized in Figure 58 through 60 for three locations 
on the MDC orbiter centerline. The most optimistic correlation would predict 
transition at 1,136 seconds and the most conservative at 384 seconds in the 
reentry trajectory for X/L = 0.5. The corresponding altitudes at which transition 
would be expected are 78 km (255,000 ft) and 57 km (186,000 ft) and were derived 
from the curves of Figures 61 through 63. These curves contain a superposition 
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FIGURE 55 

GD/C BOOSTER UNCERTAINTIES 
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FIGURE 58 


UNCERTAINTY IN TIME OF TRANSITION ONSET 
BASED ON LOCAL REYNOLDS NUMBER CORRELATION 
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FIGURE 59 

UNCERTAINTY IN TIME OF TRANSITION ONSET 
BASED ON MOMENTUM THICKNESS REYNOLDS NUMBER CORRELATION 
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FIGURE 60 

UNCERTAINTY IN TIME OF TRANSITION ONSET 
BASED ON CORRELATION OF THE MDAC PARAMETER 
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2.88 

279 

2.63 
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2.55 
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of the trajectory variations of the transition parameters, the data correlations 
and lines corresponding to one through three standard deviations of the 
transition data about these fits. Transition is predicted to occur at the inter- 
section of the trajectory and correlation curves. 

A summary of transition altitudes for the better correlation attempts described 
in Section IV are contained in Figure 64. Altitudes corresponding to nominal and 
one standard deviation are listed for X/L = 0.5, Based on these results, the 
correlations which include the ratio of wall to local temperature and surface 
deflection angle have the smallest uncertainty in altitude. The nominal transition 
altitude is between 70 km (230,000 ft ) and 72 km (235,000 ft ) for these correl- 
ations. The altitude uncertainty corresponding to one standard deviation is less 
than 4.6 km (15,000 ft ). The altitude uncertainty for the more commonly used 
correlation forms is approximately twice as great. Care should be used in general- 
izing these results since it is well known that extrapolation of ground test 
transition data to flight is still questionable. Therefore, it is recommended that 
uncertainties based on both optimistic and conservative transition be evaluated 
and compared. 

le to Flight - As previously discussed, the question of scaling 
uncertainties measured in ground test to the flight situation is not completely 
resolved. It is recognized that the ground test data contain sources of error 
which are not present in flight. Accuracy analyses by tunnel operators usually 
indicate that these error sources are small (although there is some indication 
that data taken by the Shuttle contractors using Stycast models may contain 
significant errors). If the instrumentation -related errors are small compared 
with the total data scatter, it is likely that the major source of the scatter is 
the inability of theory to match the actual aerothermodynamic phenomena. The 
uncertainty thus produced may therefore propagate to the flight situation. This 
study was hampered by the limited quantity of data on similar configurations 
spanning large Mach number/Reynolds number ranges with wall cooling (e.g.T /T ) 
comparable with flight. As a result, conclusive proof that the wind tunnel derived 
uncertainties could be extrapolated to flight could not be derived. Therefore, 
the various approaches to deriving uncertainty factors were compared for various 
confidence levels (both for one-and two-sided normal distributions) of the data 
about representative theoretical approaches. An example which corresponds to peak 
turbulent heating for the orbiter centerline is shown in Figure 65. Nominal and 
off-nominal heating corresponding to one, two, and three standard deviations of 
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the data about each theory were made for the Spalding-Chi , Beckwith and Gallagher 
and PV turbulent theories. It is important to note that differences exist even 
for the nominal case since the theories extrapolate differently to the flight 
conditions. This occurs even though the standard deviation of the data about 
best fits for each theory were essentially the same at the flight angle of attack. 
Comparison of the data with theory could not , therefore , be used to select the 
theory which best modeled the test data. 

Included on the curve are the heating rate excursions related to various 
statistical approaches and confidence levels. These include 95-percent con- 
fidence limits for the mean (indicated in the figure as DF) and confidence limits 
with tolerance for two sided confidence levels for the data sample of 166 points. 
Because of the large sample size, the factor placed on the mean would be quite 
small and, if used as a design factor, would not even encompass the range of 
heating rates predicted in flight for nominal conditions by the three theories. 

In applying tolerance limits with confidence it is recognized that the standard 
deviation is not precisely known, and, therefore, it is necessary to express 
for a given confidence level the percentage of the data which will fall in the 
interval. The intervals for 95-percent confidence that 95 and 99-percent of the 
data will fall within the interval are also included in the figure. If the 
sample were Infinite (the total population) the 95/95 location would be at 1.96 a. 
Also included on the figure are the values for 95/95 corresponding to one-sided 
confidence limits; that is, values corresponding to a 95-percent confidence 
level that the peak values will not be exceeded. 

It is interesting to note that the use of two standard deviations to represent 
a 95“percent confidence level, the approach used in the BGRV design and also 
applied by Scottaline in Reference 19 is midway between the 95/95 confidence 
limit with tolerance for one-sided and two-sided distributions. 

In Figures 66 through 68, flight test data are compared with ground test 
results as a function of angle of attack. These results for a slender axisym- 
metric configuration indicate that the three-dimensional effects are well modeled 
by the laminar ground test data (Figure 66 and 67)and that the majority of the 
flight data are within two standard deviations of the ground test scatter. The 
turbulent data (Figure 68) are slightly lower than predicted by the ground test 
derived heating method (modified Spalding-Chi) and all measurements were lower 
than predictions based on two standard deviations of the ground test results. 
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Because of limitations in the sample size and the large scatter present in the 
existing shuttle data, it does not appear justified to utilize highly sophisti- 
cated statistics to define confidence levels with tolerance since the difference 
in results from using the standard deviation as a measure of confidence is small. 
However, it is recommended that future evaluations of aerodynamic heating 
uncertainties for application to Phase C/D use the more sophisticated approach, 
one-sided tolerance limits with confidence. 
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VI THERMAL PROTECTION SYSTEM REQUIREMENTS 

Thermal protection system weights were made for the Shuttle booster and 
orbiter using the data correlations derived in this study for nominal heating and 
for uncertainties corresponding to one, two, and three standard deviations of the 
data around the nominal value. The configurations consisted of the GD/C B-15B-2 
booster and the MDC 050 B orbiter. These were selected because of the range 
and quality of the heating data available on these configurations. The trajectory 
for the booster consisted of a representative heat sink booster trajectory with a 
staging velocity of approximately 2.155 km/sec (7050 ft/sec). The orbiter 
trajectory was the MDC Phase B design baseline reentry trajectory for 2040 km 
(1100 rani) lateral cross range. The TPS concept for the booster was aluminum 
heat sink, and, for the orbiter, was RSI attached through strain isolation foam 
to aluminum structure. 

TPS requirements for the canards, control surfaces, and base region were not 
included in this study. Also, the effect of plume impingement heating to the 
booster upper surface by the orbiter engine was neglected. 

Calculations were made for all laminar and all turbulent flow, as well as for 
both optimistic and conservative boundary layer transition criteria. 

1 . Booster TPS Requirements - A sketch of the booster and the locations for 
which TPS requirements were calculated is shown in Figure 69. The vehicle length 
was assumed to be 82 meters (269 feet). The surface areas corresponding to each 
of the fourteen regions for which the calculations were made are also included 
in the figure. It was assumed that the unit weight of the TPS was uniform for 
each region. The TPS concept for the booster was aluminum heat sink, that is, 
sufficient skin thickness was provided to limit the peak temperatures of the 
aluminum structure to 450°K (350°F). The LH 2 fuel and LOX tanks were assumed 
to be an integral part of the structure , and it was assumed that internal insula- 
tion was used to maintain preflight temperatures of the LH 2 tank above the LOX 
temperature 89°K (-300°F) . This temperature was assumed to be the liftoff 
temperature in all tank areas. In the nontank areas the liftoff temperature was 
assumed to be 300°K (80°F). The thermal analysis model assumed one-dimensional 
heat flow through the skin with an adiabatic back surface. The analysis technique 
consisted of calculations of temperature histories of the aluminum for various 
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skin thickness with interpolation of the results to determine the thickness which 
would achieve a peak temperature of 450°K (350°F). 

The reference heat pulse for the booster is presented in Figure 70. The 
convective heating rates were computed for the stagnation point of a ?ne foot 
diameter sphere using the method of Fay-Riddell for equilibrium air. A wall 
temperature of 311°K (100°F) was assumed for these calculations. The peak ref- 
erence heating rate prior to separation of the booster and orbiter is approxi- 
mately 4 x 10^ watts/m^ (3.5 BTU/ft^ sec) and during entry 8 x 10 watts/m 
2 

(7.0 BTU/ft sec). The small total heat load and the fact that the booster tank 
structure has a large potential for heat absorption since it is initially at cryo- 
genic temperatures make the use of heat sink TPS a logical choice for the booster. 
A skin thickness of 0.46 cm (0.18 in.) would be able to absorb this heat load 
to a 1-ft diameter nose. This is less than the structural requirements over much 
of the vehicle. 

a. Heat Sink Requirements - The initial analyses were conducted to establish 
for the booster skin thicknesses to absorb the convective heat load independent of 
the structural requirements. The nominal heating rates were based on the correla- 
tions for the booster heating described in Section IV at the 14 locations 

shown in Figure 69. The total surface area for which thicknesses were deter- 
2 2 

mined was 3496 m (37,607 ft ). Figures 71 and 72 contain summarized unit 

weights and total heat sink weight in SI and engineering units respectively. For 

2 

the nominal heat pulses the average unit weight of heat sink is 13.6 kg/m (2.8 
2 

lb/ft ) which is equivalent to an average skin thickness of 0.51 cm (0.20 in.). 

It should be noted that this is an average, just as the tabulated values are the 
average thicknesses in the region for which the calculations were made. These 
results correspond to the minimum weight heat sink TPS, since they do not account 
for structural requirements in the tank areas, which encounter critical loads 
during boost when temperatures are low. 

Heat sink weights are also summarized for heating pulses which are 1, 2 and 
3a higher than nominal. These were obtained by applying the heating multipliers 
summarized in Figure 55 of Section V. As shown in Section V, these uncertainties 
correspond (approximately) to 68, 95 and 99.7-percent confidence levels. Using 
more refined statistical approaches would only slightly change the ra her large 
uncertainties shown in the table. Therefore, more elaborate analysis techniques 
are not justified. 
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If thermal considerations alone are looked at, the uncertainty in heat sink 
weight is approximately 85 percent for a confidence level of 99,7 percent* 

A review of the results shows that the major weight increment is on the wing lower 
surface. This occurs because of the large area involved. The average upper 
surface unit weights are considerably less than those on the lower surf ace, and thus 
have a smaller impact on the heat sink requirements. 

The above figures were based on boundary layer transition onset predictions 
using the MDAC-E developed transition criterion which was recommended by the 
Shuttle Aero thermodynamic Working Group as an interim criterion. Total vehicle 
heat sink weights were also made to assess the effect of heating uncertainties for 
assumed completely laminar and completely turbulent flow. Weight increments (from 
nominal) which resulted from these analyses are shown in Figure 73. These 
show that a weight reduction of 21,791 kg (48,041 lb) is possible if the flow 
were laminar and a weight increase of 4270 kg (9413 lb) if the vehicle were 
designed for turbulent flow. All weight increments are relative to the nominal 
weight based on the MDC Phase B transition criterion which is 45,047 kg 
(99,179 lb). 

The effect of transition criterion on the TPS weight is shown in Figure 74. 

It is interesting to note that use of the optimistic fit to the Shuttle transition 
data, i.e., local Reynolds number versus local Mach number, has only a slight effect 
on the booster heat sink weight. The transition criterion used in this comparison 
is shown in Figures 61 and 63 of Section V. The uncertainties from these figures 
corresponding to the transition data scatter for 1, 2 ,and 3a were used to evaluate 
the effect of transition uncertainty on the TPS weight. These results were combined 
statistically with the effect of uncertainty in heating rate. The combined 
uncertainties are summarized in Figure 75 from both the MDC and optimistic trans- 
ition criteria. The combined weight uncertainty was obtained as the root sum 
square of the uncertainty due to heating rate prediction and boundary layer 
transition data scatter. That is, 

AW = ^|(AW q ) 2 + (AW t ) 2 

where AW^ was obtained from Figure 74 and AW^, from Figures 71 and 72 at corres- 
ponding values of o. As shown in Figure 75, 34,981 kg 

(77,118 lb) would be required to assure a confidence level of approximately 99.7 
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FIGURE 73 

TPS WEIGHT INCREMENTS DUE TO ALL LAMINAR AND ALL TURBULENT FLOW 
GD/C BOOSTER LOWER SURFACE 


AW IN KILOGRAMS 


FLOW 

NOMINAL 

la 

2 a 

3a 
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FIGURE 74 

TPS WEIGHT INCREMENTS DUE TO UNCERTAINTY IN TRANSITION 
GD/C BOOSTER LOWER SURFACE 
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FIGURE 75 

TPS WEIGHT INCREMENTS DUE TO COMBINED UNCERTAINTIES IN HEATING 
RATE AND TRANSITION ONSET-GD/C BOOSTER LOWER SURFACE 
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percent if the MD.C transition criterion is used and 24,379 kg (53,746 lb) if a 
more optimistic transition criterion is utilized. 

b. Combined Heat Sink and Structural Requirements - It should be recognized 
that the previous discussion concerned heat sink requirements only. Because the 
booster TPS is also structure, the values In Figures 71 and 72 must be adjusted to 
reflect structural requirements. This is done in Figures 76 and 77 for SI and 
engineering units, respectively . 

Structural requirements for the LOX tank, internal tank region, and LH^ tank 
were developed in Reference 20. Using the values reported in the above reference 
it was found that the LOX tank, internal tank, and LH^ tank upper surfaces were 
sized by structural requirements. In addition, heat sink requirements were less 
than structural requirements for the nominal and la heating on the lower surface 
of the LOX and LH^ tanks. A portion of the wing upper surface was also sized by 
structural rather than heat sink requirements. Integrating these effects result 
in a considerably lower TPS uncertainty (in percent of dry weight) than if thermal 
considerations alone are included. The 84-percent uncertainty for a confidence 
level of 0.997 is reduced to approximately 55-percent of dry weight. 
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FIGURE 76 

COMBINED STRUCTURAL AND HEAT SINK WEIGHT BREAKDOWN 
NOMINAL MDAC-E TRANSITION CRITERION 
GD/C HEAT SINK BOOSTER 
S. I. UNITS 
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FIGURE 77 

COMBINED STRUCTURAL AND HEAT SINK WEIGHT BREAKDOWN 
NOMINAL MDAC-E TRANSITION CRITERION 
GD/C HEAT SINK BOOSTER 
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2 * Orbiter T PS Requirements - A sketch of the orbiter configuration used to 
define the TPS requirements is shown in Figure 78. This configuration was 
selected since it was a configuration for which a large amount of test data were 
available and represented a viable design for the fully reusable shuttle system. 

The MDC Phase B studies had baselined a metallic TPS for this configuration. More 

recent NASA decisions have led to selection of a reusable surface insulation (RSI) 
as the shuttle baseline TPS. Since the TPS weight sensitivities would be different 
for a metallic than for an RSI system, an RSI thermal protection system was 
studied to allow more direct application of these results to the current Shuttle 
orbiter. A sketch of the TPS installation is shown in Figure 79. Mullite RSI 
is attached through a foam layer to the aluminum structure. The use of mullite 
provides reuse capability to 1530°K (2300°F) and overshoot capability to 
approximately 1860°K (2900°F). The foam layer is applied to isolate the RSI from 
the structural bending loads carried by the aluminum structure,, thereby pre- 
venting potential damage to the RSI by structural buckling. The sizing of the RSI 
thickness and foam had two temperature constraints: an RSI-to-foam-bondline temper- 
ature limit of 590°K (600°F) and a peak aluminum structural temperature of 450°K 
(350 F) . Because of handling and fabrication limits the minimum thickness of RSI 
which can be applied is 0.635 cm (0.25 in.). The weights presented herein include 
the weight of RSI, its silicone waterproof coating, strain isolation foam, and 
adhesives required to bond the elements of the thermal protection system together. 
An aluminum thickness of 0.15 cm (0.06 in.) was assumed in these analyses. 

The RSI was sized for reentry only, since the small boost phase heat pulse will 
not produce critical bondline or structural temperatures and this heat will be 
dissipated during orbit for the majority of Shuttle orbiter missions. The possible 
exceptions to this case are abort or once-around missions. However, even for these 
remote situations, the boost phase is not likely to have a major effect on the 
sensitivities derived from this study for the reentry heat pulse. The initial 
temperature of the TPS component was assumed to be 311°K (100°F) , a value consistent 
with normal vehicle attitude and orbital thermal control and orbit inclination. 

Typical reentry heat pulses are included in Figure 80 for a lower surface 
centerline location at 10 percent of the vehicle length. The three heat pulses 
correspond to all laminar, to all turbulent, and to assumed transition onset based 
on the MDC Phase B transition criterion (Figure 47, Section IV). For the 
transitional heat pulse, fully developed turbulent flow was assumed to occur when 
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FIGURE 80 

REENTRY HEAT PULSE 

MDAC DELTA WING ORBITER LOWER SURFACE CENTERLINE 
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transition onset was predicted at X/L = 0.05. That is, the extent of transitional 
flow was equal to the length of the vehicle for which the flow was laminar. The 
reentry time from 122 km (400,000 ft) is approximately 1500 seconds for this 2040 km 
(1100 nmi) cross range trajectory. The heating rate histories were computed using 
the correlations developed for the orbiter in Section IV. For the lower surface 
centerline these consisted of heating rates computed using the Eckert Flat Plate 
Reference Enthalpy Method with local flow conditions based on equivalent cone flow. 
The analyses were made for real air in thermodynamic equilibrium. Correlation 
factors. Figure 34 of Section IV, were applied to these calculations as a function 
of angle of attack to account for the three dimensional characteristics of the 
flow. For turbulent flow a similar approach was utilized. The heating was computed 
USin -2/3 G Spalding-Chi corre lation with a Reynolds analogy factor corresponding to 
(P r > ' Factors developed from the test data correlations to account for 
three-dimensional flow field effects were also applied to these heating rates. 

The correlation curve used for the lower surface centerline is also included in 
Figure 34 of Section IV. 

TPS requirements were determined from correlations of TPS unit weight as a 
function of integrated heat load. The average unit weight and total TPS weight 
are summarized in Figures 81 and 82 for SI and engineering units, respectively . 

The weight breakdown is summarized for the eight sections into which the vehicle 
was divided. The nominal heat pulses and heat pulses corresponding to heating 
multipliers derived from Figure 56 of Section V for 1,2, and 3a uncertainties 
were used to define these TPS requirements. The MDC transition criteria was used 
m these analyses. The results show that the orbiter TPS weight is much less 
sensitive to the aerodynamic heating uncertainties than booster TPS weights due 
to the fact that the majority of the aerodynamic heating is reradiated to space. 

As the heating rate is increased, a greater fraction of the total heat is 
reradiated. Thus, the percent which is conducted to the structure is reduced. 

A 3a design would require a 15-percent weight increase. The greatest weight penalty 
associated with the heating uncertainties is on the wing lower surfaces, 
primarily because of the large area involved. • The location, having the largest 
uncertainties in heating are regions for which the nominal heating is quite low; 

it was found that only slight increases above minimum RSI thickness were required 
for the 3a case. 

Figure 83 summarizes the weight increments (from the nominal 
in Figures 81 and 82) if laminar or turbulent flow is assumed for the entire 
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FIGURE 83 

TPS WEIGHT INCREMENTS DUE TO 
ALL LAMINAR AND ALL TURBULENT FLOW 
MDAC DWO LOWER SURFACE 
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flight. The small increment for the nominal turbulent case results from the fact 
that the lower surface of the orbiter is essentially designed by fully turbulent 
flow if the MDC Phase B transition criterion is used. This figure shows that 
a significant weight reduction is possible if the flow remains laminar in flight. 
Assessment of the effect of the transition criterion and its uncertainty was made 
and is summarized in Figure 84. This figure shows that the use of an optimistic 
transition criterion can reduce nominal TPS requirement by approximately 2240 kg 
(5000 lb). However, if the effect of scatter of the data about the best fit 
to the transition data is includedjthe difference in TPS weight becomes insignifi- 
cant. This is shown graphically in Figure 85, in which the unit weight of the 
TPS is shown for the best correlation fits and uncertainties associated with 1, 2, 
and 3a of the transition data about the best fits. At this location, which 

corresponds to an X/L of 0.5, the maximum difference (conservative to optimistic) 

2 ' 2 

in unit weight is reduced from 3.4 kg/m (0.7 lb/ft ) for the nominal values to 
2 2 

0,39 kg/m (0.08 lb/ft ) for a 3a deviation in transition onset. 

The combined effect of heating prediction and transition onset uncertainties 
are compared in Figure 86 for the MDC and optimistic transition criteria. These 
were combined using the root sum square technique discussed above. This figure 
shows that the 3a design would be 2050 kg (4510 lb) heavier than the nominal using 
the MDC criterion and that the 3a design would actually be 91 kg (200 lb) 
lighter than the nominal based on the MDC transition criterion. This shows that 
the selection of the transition criterion to apply to the orbiter has a major 
influence on TPS weight, a situation which was not critical for the booster. The 
use of the optimistic criterion would result in a 3a design approximately 13 
percent lighter than would be predicted using the MDC criterion. 

The above weight studies and tradeoffs were made under the implicit assumption 
that la corresponded to a 68-percent confidence limit that the heating would not be 
exceeded in flight. As discussed in Section V, this use of the standard deviation 
to define confidence limits is not precise. Therefore, the more sophisticated 
approach which develops confidence limits with tolerance has been used to estimate 
the effects of heating uncertainties on orbiter TPS weight. These results 
are summarized in the detailed weight breakdown of Figures 87 and 88 for SI 
and engineering units respectively. To achieve a 95— percent confidence level that 
95— percent of the data in flight will not exceed the design requirements (95/95) a 
13“P ercen t increase in TPS weight is required. By comparison with the results in 
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FIGURE 84 

TPS WEIGHT INCREMENTS DUE TO UNCERTAINTY IN TRANSITION 
MDAC DWO LOWER SURFACE 
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FIGURE 85 

EFFECT OF TRANSITION UNCERTAINTY 
ON TPS UNIT WEIGHT 
MDAC DELTA WING ORB I TER 
X/L = 0.5 
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FIGURE 86 

TPS WEIGHT INCREMENTS DUE TO COMBINED 
UNCERTAINTIES IN HEATING RATE AND TRANSITION ONSET 
MDAC DWO LOWER SURFACE 
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Figure 81, a corresponding 2(7 design would require a 10-percent increase in 
weight. Thus, the use of more rigorous statistics results in slightly higher heat- 
ing multipliers than derived from direct use of values at 1,2, and 3cr . 

Since the upper surface of the wing encounters low heating and the TPS is 
sized entirely by the minimum RSI thickness, a study was made to evaluate the use 
of thickened aluminum or titanium heat sinks. Unit weights of these heat sink 
concepts are compared with the RSI TPS as a function of trajectory heat load in 
Figure 89. This curve shows that a 1<T design is possible with no weight penalty 
for the aluminum heat sink concept, and a 2cr design if titanium is used. However, 
the RSI concept provides the lightest weight if it is necessary to design for a 
3(7 heat pulse. 
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VII CONCLUSIONS AND RECOMMENDATIONS 

This study has provided an assessment of the effect of aerodynamic heating 
uncertainties on Space Shuttle TPS weights, . thereby identifying the prediction 
methods for which improvement is most important. These include windward turbulent 
heating for the booster and orbiter and development of more reliable boundary layer 
transition criteria for the orbiter. Mated interference heating was found to be 
less critical than initially suspected because of the low staging velocity of the 
heat sink booster and apparent Mach number dependence of the interference heating. 
Uncertainties in the heating to upper surface shielded regions at reentry angles 
of attack were quite large. However, for the vehicle and TPS concepts studied in 
this program, the heat sink or RSI thickness in shielded regions was essentially 
sized by structural or manufacturing limitations. These uncertainties could, on 
the other hand, be critical for other configurations or TPS concepts. 

The heat sink booster was found to be much more sensitive to aerodynamic 
heating uncertainties and a high confidence (3a) design would be approximately 
33-percent heavier (in skin and structure weight) than required if nominal 
heating is assumed. The equivalent penalty in TPS for the orbiter would be 
approximately 13-percent. In addition to the above general conclusions, a 
number of specific observations can be made: 

1, The Shuttle heating data base is reasonable for the current 

state of Shuttle development. However, this study was hampered by certain limita- 
tions in this data base. For example, it was difficult to separate the effect of 
geometric and environmental variables on the heating uncertainties , since few 
configurations were tested over a wide range of Mach and Reynolds numbers. In 
particular, the variations in tunnel conditions prevented the selection of a 
theory which best matched the data even though differences in scale to flight would 
result for each theory. In addition, the emphasis on paint data made the data 
reduction and evaluation process difficult and increased the derived uncertainties. 
Care was taken to minimize these effects, but the reported uncertainties do 
contain some. 

2. Mated interference heating uncertainties, though large, were not a major 
contributor to the weight uncertainties because of the low staging velocities and 
correspondingly low interference heating multipliers predicted by the correlation 
used in this study. Alternate correlation approaches, in terms of Reynolds number 
or both Mach number and Reynolds number, could affect this conclusion. These 
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were not tried because of the limited amount of low Mach number data on the mated 
configurations . 

3. The statistical approach to evaluating data scatter and applying the 
results to flight were reviewed. A more sophisticated approach (than simply 
relying on the standard deviation of the data sample) to establish confidence 
levels was developed. This approach provides confidence that a percentage of 
the data population will be less than the design value with some confidence level. 
This approach is valid if it is assumed that the factors affecting ground test 
scatter will also occur and are scalable to flight. It is known that this is not 
a universal truth , since a portion of the ground test scatter may be due to data 
fits or instrumentation error. However, if the major portion of the scatter is due 
to improper modeling of the data by theory, such scaling is possible. Unfortun- 
ately, because of insufficient flight test data, statistics cannot be employed 
to validate whether ground test data and their uncertainties scale to flight. 
However, comparisons of the scatter in BGRV flight data with ground test data 
and predictions provide some assurance that the approach is reasonable. 

4. The uncertainty in predicting lee side heating using simple correlating 
approaches was found to be quite large at the reentry angles of attack. Simple 
correlation m terms of h/h^ produced a scatter of approximately an order of 
magnitude. More careful screening of the data or improved correlation should 
reduce this uncertainty. Further development of lee side methodology is certainly 
needed, even though these uncertainties were not found to be critical for the 
fully reusable Shuttle concept. Should it be desired to fly at lower angles of 
attack during reentry, the heating level would increase and large uncertainties 
could penalize the design. For that matter, if the orbiter TPS concept were 
changed to heat sink or heat shield on shielded areas, large weight penalties 
would be required to design for these uncertainties. It should be noted that the 
low tunnel heating rates measured in these areas make the data sensitive to 
instrumentation errors. More data are necessary to reduce the shielded region 
uncertainties, but it may be necessary to improve tunnel instrumentation if 
reliable data are to be obtained. 

5. The boundary layer transition criterion was found to be a significant 
contributor to the TPS weight uncertainty for the orbiter. This uncertainty will 
probably not be greatly reduced through additional ground tests. However, 
additional ground testing can provide an understanding of the sensitivity of 


128 


MCOO^\fLL DOUGLAS 4STffOA/4L/r/CS COMP4/VV - C 49 T 



30 JUNE 1972 


MDC E0639 


transition to configuration variables, e.g., nose geometry, lower surface curvature. 
It should also be noted that the better correlating approaches used a configuration 
related parameter, Re@* The cooling parameter which was found to improve the 
correlation was the ratio of wall to edge temperature rather than wall to recovery 
temperature. For this correlation, transition on the orbiter was predicted at 
72 km (235,000 ft), whereas, the MDC Phase B criterion predicted transition at 
78 km (255,000 ft) and correlation of the data in terms of local Reynolds number 
at 58.5 km (192,500 ft). 

6. Except at angles of attack near zero, the use of normal shock entropy 

to predict local flow conditions provided no improvement in the data correlation. 

The edge properties used in this study were analytically derived for cones having 
equivalent deflection angles to the local surface angles. These were computed only 
for those cases where the data were compared with a theoretical two-dimensional 
prediction method. For data correlated to Fay-Riddell stagnation point theory, 
such assumptions are not necessary. Since the Shuttle configurations are quite 
complex, it is probable that the data scatter could have been reduced if the 
boundary layer edge properties had been better defined. This is particularly true 
for the turbulent data which is more sensitive to local Reynolds number than for 
the laminar case. 

7. Many data points were not evaluated in this study either for lack of 

resources or unavailability of the data in time for use. The tests have been 

cataloged and are reposing in SADSAC. It would seem desirable to develop automated 
processing techniques which would allow correlation of these results without the 
time consuming reading of isotherms or hand conversion of coefficients to appro- 
priate equivalent theoretical local flow conditions. Also, the data in SADSAC 
should be screened for errors, and invalid or inaccurate data removed. 

8. It is recommended that particular care be taken in Phase C/D that (1) a 

statistically meaningful data base be obtained, (2) instrumentation errors be 

assessed, and (3) best methods for scaling to flight be developed on the basis 

of the match of the test data to theory. 
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VIII LIST OF SYMBOLS 


Symbol 

C 

H(T ), H(.9TO), h 
K 

L 

m 

M 

N 

n 

P 

q 

Re 

Re^/L 

Re e 

S 

5 

SI 

T 

V 
w 
X 
x 

Y 

a 

6 

9 

R 


Definition 

Model Wing Local Chord Length 
Heat Transfer Coefficient 

Multiple of the standard deviation corresponding to various 
proportions included within the tolerance range 
Model reference length (body axial length) 

Meter 

Mach number 
Normal distribution 
Sample size 
Pressure 
Heating rate 
Reynolds Number 

Free-stream Unit Reynolds Number 

Momentum Thickness Reynolds Number 

Sample standard deviation 

Model Wing Exposed Semi-span Length 

Standard International 

Temperature 

Velocity 

TPS weight 

Axial distance 

Sample mean 

Wing outboard distance from root 
Angle of attack 
Local Flow Deflection Angle 
Shock angle 
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Symbol 

A 

U 

* 

X 2 

a 

a 

t 

A 

Subscripts 

B 

B+O 

£ 

e 

L 

tr 

♦ 

au 

s 

Ref 


w 


Definition 

Leading edge sweep angle 
Population mean 
Peripheral angle 
Chi-square distribution 
Population standard deviation 

Unbiased estimate of the population standard deviation 
Time, student’s t distribution 
Increment in parameter prefixed 

Booster alone 

Booster mated with orbiter 
Lower surface centerline value 
Edge conditions 
Local conditions 
Transition location 
Local peripheral value 
Free stream 
Stagnation point 

Reference value of h or q for a unit (R=l f ) sphere 
scaled to model size. 

Wall condition 
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